Fermionizing a small gas of ultracold bosons 
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We study the physics of a rapidly-rotating gas of ultracold atomic bosons, with an internal degree 
of freedom. We show that in the limit of rapid rotation of the trap the problem exactly maps onto 
that of non-interacting fermions with spin in the lowest Landau level. The spectrum of the real 
bosonic system is identical to the one of the effective fermions, with the same eigenvalues and 
the same density of states. When the ratio of the number of atoms to the spin degeneracy is an 
integer number, the ground state for the effective fermions is an integer quantum Hall state. The 
corresponding bosonic state is a fractional quantum Hall liquid whose filling factor ranges in the 
sequence i/ = 1/2,2/3,3/4,..., as the spin degeneracy increases. Anyons with 1/2,1/3,1/4,... 
statistics can be created by inserting lasers with the appropriate polarizations. A special situation 
appears when the spin degeneracy equals the number of atoms in the gas. The ground state is 
then the product of a completely antisymmetric spin state and a = 1 Laughlin state. In this case 
the system exhibits fermionic excitations with fermionic statistics though the real components are 
bosonic atoms. 



I. INTRODUCTION 

Ultra cold atomic gases have been shovifn to be ideal 
candidates to observe several intrinsically quantum phe- 
nomena, since they can be easily controlled and manip- 
ulated by electromagnetic fields. What type of quantum 
phenomena we can observe depends dramatically on the 
strength of the atomic interaction, going from one parti- 
cle phenomena to highly entangled quantum states. 

The first experimental highlight with these systems 
was the achievement of Bose-Einstein condensation of 
weakly interacting atomic gases [Q. This success mo- 
tivated a lot of experimental and theoretical work with 
cold bosons in the limit of weak interactions. In this 
regime atoms essentially occupy the same single particle 
state, and the physics of the system is perfectly under- 
stood via a Gross-Pitaievski equation for the macroscopic 
wave function. Many interesting phenomena have been 
studied in this regime and there still exist certain open 
questions ||^. 

On the other hand, it seems clear that going in the 
direction of strongly correlated atoms, and exploring 
what happens beyond the validity of mean field theory, 
promises to open a new path to the observation of novel 
physical phenomena. Strong correlations between parti- 
cles have been shown to lead to unusual fascinating ef- 
fects in condensed matter systems. In the same way, if 
an atomic system is driven into the strongly interacting 
regime, similar interesting phenomena may be observed. 
Moreover, since atomic gases can be well controlled and 
manipulated, they may provide us with new scenarios for 
the observation of these phenomena and, perhaps, new 
ones that have not been accessible so far in solid state en- 
vironments. The difficulty here is, however, how to drive 
an atomic system into a strongly interacting regime. In 
order to enter this regime one should make the typical 
interaction energy larger than any other (single particle) 
energy scale in the problem. Naively speaking, it may 



seem that the easiest way to do this would be just to in- 
crease the interaction strength, either by increasing the 
scattering length or the atomic density. However, though 
the scattering length can be effectively enlarged by using 
Feshbach resonances and the atoms squeezed to in- 
crease the density, both methods may lead to undesirable 
instabilities in the system. 

Alternatively, in order to make interactions to dom- 
inate, one can decrease the single particle energies by 
creating degeneracies. Here, one possibility is to use in- 
ternal levels, as in the recent proposal to entangle atomic 
beams by generation of spin-squeezed states Q . Degen- 
eracies can also be induced by using a periodic external 
potential, since single particle states corresponding to dif- 
ferent wells will have the same energy. In the limit in 
which the amplitude of tunneling between wells is small 
compared to the typical interaction energy, the system 
has been predicted to undergo a transition from a su- 
perfluid phase to a Mott insulator phase which is 
highly correlated. This quantum phase transition has 
been recently achieved experimentally in an optical lat- 
tice . Another way of achieving a quasi-degeneracy in 
the atomic motional states is to rotate the trap that con- 
fines the atoms. When the frequency of rotation is small 
compared to the frequency of the trap, the system forms 
a vortex condensate, with a number of vortices that in- 
creases as the rotation becomes faster ||^. This vortex 
condensate is a non-correlated state, in which all the 
atoms occupy the same single particle quantum state. 
However, when the rotation frequency is close to the 
harmonic trap frequency the system has been predicted 
to enter a fractional quantum Hall regime, in which the 
atoms are strongly correlated ^. The ground state is 
a ^-Laughlin state, a highly entangled liquid with nearly 
uniform density. This atomic liquid offers the fascinat- 
ing possibility of creating and manipulating anyons with 
i-fractional-statistics 0. By piercing the system with 
off resonant localized lasers, i-Laughlin quasiholes may 
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be created and manipulated to probe directly their frac- 
tional statistics. For the fractional quantum Hall regime 
to be experimentally achieved in current rotating traps 
one would need bosonic gases consisting of a small num- 
ber of atoms, a limit that has remained practically un- 
explored. A particularly interesting exception is a recent 
experiment with an optical lattice that has achieved 
a regime in which a few number of atoms are confined in 
each well of the lattice. Since the two-dimensional rotat- 
ing realization of this lattice may be accessible in near 
future experiments, further theoretical investigation in 
the physics of small rotating bosonic gases is highly mo- 
tivated. 

In this paper we extend the results of our previous 
work 1^ on a rapidly rotating gas of bosonic atoms, to 
the case in which the atoms have an internal degree of 
freedom. We will show that the addition of internal levels 
introduces new aspects of the problem, and novel exotic 
correlated phases come out. 

The problem of strongly interacting bosons with spin 
in a rapidly rotating trap is, in principle, very difficult 
to tackle; in particular it is not easy to guess what type 
of ground states and excitations appear. In this work we 
present a formalism, what we call fermionization of the 
bosonic system, within which the problem takes a very 
simple form. We will show that in the limit of rapid rota- 
tion and small temperatures, the problem of interacting 
bosons maps exactly onto a problem of non-interacting 
fermions. The ground state corresponds to a Fermi sea 
in which the fermions fill up the single particle states 
in order of increasing energy. Furthermore, the lowest 
energy excitations consist of fermion-hole excitations in 
the vicinity of the Fermi level. The spectrum of the real 
bosonic system is identical to the one of the effective non- 
interacting fermions, with the same eigenvalues and the 
same density of states. 

The fermionization scheme is valid for any number of 
particles N and any number of internal states n, pro- 
vided that the atomic interaction does not depend on 
the internal states of the atoms. Perhaps the most inter- 
esting states appear whenever the ratio N/n is an inte- 
ger number. The effective fermions then form an integer 
quantum Hall state with filling factor ly = n [|lOj. The 
corresponding bosonic state is a fractional quantum Hall 
liquid with filling factor v — n/(l + n). This fractional 
liquid is a multicomponent state made out of n copies of 
a Y^~Laughlin-liquid one for each internal state. 
Thus, by increasing the number of internal states we can 
go beyond v = 1/2, and lead the system into new frac- 
tions in the sequence — 1/2, 2/3, 3/4, .... All type of 
anyons with statistics p/q, p < q can now be created 
by choosing the adequate spin degeneracy and inserting 
lasers with a polarization and frequencies which couple 
to a given atomic internal state (spin component). 

An exotic situation appears when the number of inter- 
nal states equals the number of atoms in the trap. The 
effective fermions form in this case a pure spinor state 
with no spatial dependence. The corresponding bosonic 



state is then the product of a completely antisymmetric 
spin state and a, v = 1 Laughlin state, which is fully an- 
tisymmetric too. This means that as long as the spin 
degrees of freedom are not excited the behaviour of the 
system will be fermionic. We will show that the elemen- 
tary excitations of the systems are true fermions with 
fermionic statistics, though the real components of the 
system are bosons. 

The formalism of fermionization that we present in this 
work can be also applied to explore several interesting 
possibilities. A very attractive one consists of inducing 
pairing in the effective fermions, so that they form a su- 
perconducting state. We will show that this is in fact 
possible when the atomic interaction depends on the in- 
ternal state. We will consider the case of atoms with 
two internal levels when the interaction between atoms 
in the same spin state is much stronger than the one be- 
tween atoms in a different spin state. We will see that the 
effective fermions configurate in pairs, forming a paired 
state identical to the BCS state of electrons with spin 
triplet p-wave pairing. The corresponding bosonic state, 
a Pfaffian state, is very well known in the context of the 



quantum Hall effect |13 



This paper is organized as follows. In section II we de- 
velop the fermionization scheme that maps the problem 
of interacting bosons in a rapid rotating trap onto one of 
non-interacting fermions. Using this formalism we show 
in section III that, when the ratio N/n is an integer num- 
ber, the ground state of the system is a v = n/{l + n) 
fractional quantum Hall liquid. Starting with these corre- 
lated liquids, we show how to create anyons withp/(l-|-n) 
fractional statistics, with 1 < p < n. We also discuss 
a possible experiment to reveal directly the statistics of 
anyons in a Ramsey-type interferometer. In section IV we 
study the special situation in which N = n. In section V 
we illustrate the results of the preceeding sections with 
some numerical exact diagonalizations for N = 4,6,8 
particles. The experimental conditions required to ob- 
serve the fractional quantum Hall states and the anyons 
described in this work are analyzed in section VI. As an 
example of the interesting possibilities of the fermioniza- 
tion scheme we show in section VII how pairing in the 
effective fermions can be induced in a system of atoms 
with two internal levels. Finally, in section VIII we sum- 
marize the results of this work, and discuss a possible 
scenario to observe the fractional quantum Hall states 
and the anyons. 



II. THE FERMIONIZATION 

We consider a set of N bosonic atoms confined in a 
potential which rotates in the x — y plane at a frequency 
Q,. The atoms have an internal degree of freedom, that 
we will call spin, which can take n = 2s -I- 1 values. We 
will assume that the confinement in the z direction is suf- 
ficiently strong so that we can ignore the excitations in 
that direction, and we will consider that the resulting po- 
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tential in two dimensions is isotropic and harmonic. The 
Hamiltonian describing this situation in a frame rotating 
with the trap is: 

i=l ^ ' i<j 

(1) 

with Liz being the z component of the angular momen- 
tum of the i-th atom, and where we have used the trap 
energy, fkj, as the unit of energy, and i = {h/mwY^'^ 
as the unit of length. The atoms are interacting via an 
effective contact potential, with an interacting coupling 
constant g related to the s-wave scattering length, a, 
and to the localization length in the z direction, i^i by 
g — y^2/TTa/£z. In most parts of this paper we will also 
assume that the interaction between the atoms is spin 
independent. 

Following our previous work we begin by writing 
the Hamiltonian (|l|) in the following form: 

H = Hb+Hl + V. (2) 

Here, Hb = J2iLi "^f/^+J^i /2— is the quantum Hall 
single particle Hamiltonian, whose single particle energy 
levels are the Landau levels equally spaced by the cy- 
clotron energy 2hu!. The Hamiltonian Hl = (1 — fl/uj)Lz 
is proportional to the z component of the total angular 
momentum, Lz = X^il^i^iz, and V is the interaction 
term. 

From now on we consider the rapid rotation limit, in 
which the energy scales characterizing Hamiltonians Hb 
and V are much larger than the one corresponding to Hl ■ 
In this limit the ground state and elementary excitations 
of the system must lie on the subspace of common zero 
energy eigenstates of Hb and V. 

We will show that, projected to the subspace J^, the 
bosonic problem described above exactly maps onto a 
problem of non-interacting fermions with n internal lev- 
els. We will first show that any wave function in T is 
the product of a = 1 Laughlin state and a completely 
antisymmetric wave function in the lowest Landau level. 
We will then see that this antisymmetric wave function 
describes an effective system of non-interacting fermions. 

Let be a wave function in J^. In order to be an eigen- 
state of Hb with eigenvalue zero, ^E" must lie within the 
subspace spanned by tensor products of single particle 
states in the lowest Landau level In the symmetric 
gauge these single particle states are labeled by the third 
component of the angular momentum m and the spin a, 
and have the following form: 

where ^ = {z,u), z = x + iy, and is a s-spinor with 
components u'^iP) = Sg-fj, for cr, /3 = — s, ... , s. It follows 



that must have the form: 

mi ,. . . ,771^ 
CTi ,crjv 

with nii = 0, . . . , cxD, ai = — s, ... ,5, and where for 
simphcity we have omitted the global exponential factor 

The wave function ^E* must be an eigenstate of V with 
eigenvalue zero as well. We will show that the function 
ni<j('^i — Zj) is a factor of ^. Let's choose any pair of 
particles i and j. The dependence of 4" on Zi and Zj 
can be reexpressed in terms of the relative and center 
mass coordinates, z.y, Z.y, so that we can expand ^'[^] — 
Sm ^ij^rn, whcrc Fm dcpcuds on Zij, Ui, Uj, and on the 
positions and spinors of all the other particles. In order 
for ^E* to be annihilated by the hard-core interaction V , 
Fq must be identically zero. It follows that Vi,j, is a 
factor of so that 

m = ml[('^-^j), (5) 

i<j 

where $[^] is an analytic function of Zi for all i, and 
therefore it lies in the lowest Landau level. 

The factor H[z] — Y[i<ji^i ~ ^j) proportional, up to 
a normalization constant, to the v = 1 Laughlin wave 
function. Since this factor is fully antisymmetric, the 
function (f> must be completely antisymmetric too, for 
^E* to be a bosonic wave function. The factorization (|^) 
thus states that any bosonic wave function lying within 
the subspace J-^ is completely specified by a fermionic 
wave function in the lowest Landau level. 

We will derive now an effective Hamiltonian for the 
fermionic wave function <E>, in the sense that the spectrum 
of this effective Hamiltonian will be identical to the one 
of Hamiltonian (||), projected onto the subspace J-. 

Let vE'n = 'E'nH be an eigenstate of Hamiltonian (|^) 
with eigenvalue This state satisfies 

i/*„ =. Hl^^ = (6) 

Since H is an eigenfunction of H^ with eigenvalue (1 — 
Vl/uj)N{N — l)/2, we can write the above expression in 
the form: 

eau = (-ffL$„)n + $„(i7LH) (7) 

= {l-n/uj)[{Lz + N{N -l)/2)<Pr.]Ti- (8) 

We have then proved that if \E'„ is an eigenstate of Hamil- 
tonian (H) the corresponding fermionic state <E>„ is an 
eigenstate of the effective Hamiltonian 

H,fj^{l-n/u;)[Lz+N{N-l)/2], (9) 

with the same eigenvalue. This means that our prob- 
lem of N interacting bosons projected to the subspace 
J- reduces to an effective problem of N free fermions 
in the lowest landau level, governed by the Hamiltonian 
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(^. The spectrum of the real bosonic system is iden- 
tical to the one of the effective free fermions, with the 
same eigenvalues and the same density of states. The 
fermionic eigenstates are connected to the bosonic ones 
by the transformation (^). Note that this transforma- 
tion is not unitary, and therefore an orthonormal set of 
fermionic states will be transform into a set of bosonic 
states that will not be, in general, orthonormal. 

The problem of free fermions in the lowest Landau level 
with Hamiltonian can be easily solved. The ground 
state is a Fermi sea, in which the effective fermions fill up 
the single particle states i-Pma in order of increasing angu- 
lar momentum m. Note that for each angular momentum 
m we have n single particle degenerate states. For a given 
number of particles N we can write N = kn + r, with 
< r < n. This means that we will have the first k sin- 
gle particle levels completely occupied and the fc -|- 1 level 
partially occupied with r fermions. The ground state 

will have degeneracy ^ ^ ^ ' corresponding to the differ- 
ent choices for the internal states of the last r fermions. 

If <i>i?5 is the wave function corresponding to the Fermi 
sea, the bosonic ground state is then just: 

m]^<^Fsml[(z^~z,). (10) 

i<j 

Within this effective fermionic picture it is also very 
simple to find the low-lying excitations of the system. 
They are fermion-like excitations in which an effective 
fermion crosses the Fermi level, leaving a fermionic-hole 
in the Fermi sea. 



III. MULTICOMPONENT ATOMIC 
FRACTIONAL QUANTUM HALL LIQUIDS AND 
THEIR ANYONS 

In this section we will study the case in which the 
ratio N/n is an integer number fc, so that we have a 
non-degenerate ground state with k completely occupied 
fermionic levels. We will show that the effective fermions 
form integer quantum Hall states at filling factor v = n. 
The corresponding bosonic states are v — n/{l + n) frac- 
tional quantum Hall liquids with 1/n-anyons. 

A. The u = 1/2 Laughlin liquid 

We start out with the trivial case in which n = 1 so 
that there is no spin degeneracy. The effective fermions 
will in this case occupy the single particle states with 
angular momentum m = 0, . . . , — 1. The Slater de- 
terminant corresponding to this situation has the simple 
form: 

N 

^l[z]^l[iz^^Z,). (11) 
Kj 



This state is the v = 1 Laughlin state, the same as for 
electrons in the quantum Hall effect. The density of this 
state is nearly uniform, with one fermion on average per 
unit area. The corresponding bosonic state, 

N 

nz]^X{{z^^Zjf, (12) 

i<3 

is the V — \/2 Laughlin liquid that we studied in our 
previous work Starting from this liquid one can cre- 
ate 1/2-atomic-quasihole excitations with 1/2-fractional 
statistics. The fractional statistical phase of these quasi- 
holes may be revealed directly in a Ramsey-type inter- 
ferometer 1^. 

B. The = 2/3 hquid 

Before considering the general case of any number of 
internal states, it is illuminating to study the case of two 
internal states (and an even number of atoms). We will 
use the spin language, and say that we have atoms with 
spin s = 1/2 that can be either up (t) or down (|). 

The effective fermions will fill up the single particle 
states with angular momentum m = 0, . . . , N /2 — 1 and 
both spin states a =t;i- We will have N/2 fermions 
with spin up and N/2 fermions with spin down, and each 
group will form a, v = 1 Laughlin state. The wave func- 
tion corresponding to this situation is 

,^n)}, (13) 

where 

N/2 

rA^] = n(^^ - ^ =T> I (14) 

i<j 

is a = 1 Laughlin state of N /2 fermions in the spin state 
cr, and the operator A antisymmetrizes over all possible 
ways of distributing the fermions in two groups of N/2 
fermions. The state (^3|) \s a. v = 2 quantum Hall state, 
made out of two v — 1 states, one with spin up and the 
other with spin down. The density profile of this state is 
nearly flat in the bulk, with two fermions per unit area on 
average, one with spin up and the other with spin down. 

The corresponding bosonic state is obtained from the 
wave function ( p^ multiplying by the factor Y\'^ (zi — Zj). 
To gain insight about the properties of this bosonic state 
it is very convenient to use the language of the quantum 
Hall effect, that relates angular momenta, filling factors 
and densities | p^ . In general, given an incompressible 
liquid in the lowest Landau level, its filling factor u is 
related to its total angular momentum by L ~ N"^ /2v. 
Furthermore, the liquid has a uniform density in the bulk, 
with V particles per unit area on the average. 
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Multiplication by the factor 11 increases the total an- 
gular momentum Lp oi the original fermionic state, so 
that the resulting bosonic state has a total angular mo- 
mentum 



L = Lf 



N{N -1) 



(15) 



This means that the new bosonic state will have a filling 
factor 



Up 



1 + iyp 



(16) 



As i^F = 2 we have that the bosonic state is a = 2/3 
fractional Hall liquid. This fractional liquid is made out 
of two z/ = 1/3 liquids, one with spin up and the other 
with spin down. The average density in this state is 
of 2/3 atoms per unit area, 1/3 with spin up, and 1/3 
with spin down. Roughly speaking, we can say that one 
bosonic atom is made out of three effective fermions, or 
that one fermion corresponds to 1 /3 atom. This counting 
will be very useful to study the anyonic excitations in the 
fermionic picture. 



C. The = n/{l+n) multispinor liquids 

After the above discussion the generalization to any 
number of internal states is straightforward. The ef- 
fective fermions will fill up the single particle states 
with angular momenta m = 0, . . . , N/n — 1, and spin 
a = — s, ... , s. We will have n groups of N/n fermions, 
one for each spin state, and each group will form a, v = \ 
polarized liquid of the form 



N/r, 



(17) 



The fermionic ground state is the product of the wave 
functions of the different groups, antisymmctrized over 
all possible ways of distributing N particles in groups of 
N/n particles: 



(18) 



The state ( p^ is & v = n quantum Hall state, con- 
sisting of n polarized v = 1 liquids, one for each possible 
spin state. In this state we have n fermions per unit area, 
one in each of the spin states. 

The relation ( p^ tells us that the corresponding 
bosonic state will be a fractional state with filling fac- 
tor V = ri/(l ~|- n). This liquid is made out of n polarized 
Vcr = 1/(1 + '^) fractional liquids, one for each spin state 
a. We will have n/{l-\-n) atoms per unit area, 1/(1 -I- n) 
per each spin state. 

By increasing the spin degeneracy n = 1,2,3,..., one 
can drive the atomic system into fractional quantum Hall 
states in the sequence v = 1/2, 2/3, 3/4, .... 



D. Anyons 

We will show now how to create anyons starting from 
the fractional atomic quantum liquids presented in the 
preceeding subsection. Following the ideas of our previ- 
ous work H] we insert a laser localized (within an area 
P) at some position 77. We require I77I to be much smaller 
than the size of the quantum Hall state (~ i'~^\/N — 1), 
so that the border effects are negligible. We will also as- 
sume that the laser only couples to atoms with spin ctq. 
The presence of the laser can be described by a localized 
repulsive potential, so that the new Hamiltonian for the 
bosonic atoms can be approximated by: 



(19) 



where the operator projects the spinor of the i-th 
particle on the state with spinor u'^° . In the same way as 
we did in the previous section we can project Hamiltonian 
(|l9| ) onto the subspace of wave functions of the form 
= $ Y[i<j i^i ~ ^^"i derive an effective Hamiltonian 
for the fermionic wave function $, which has now the 
form 111: 



(20) 



We will find the ground state of the fermionic Hamilto- 
nian ( po[ ) and show that the corresponding bosonic state 
is a fractional quasihole state localized at position 77. We 
consider the limit in which the intensity of the laser is 
much larger that the energy scale characterizing Hamil- 
tonian Heff. In this limit the new ground state $r)<To 
must be a zero eigenstate of the potential created by the 
laser. It follows that, whenever the z~th particle is in the 
spin state ctq , {zi — rf) must be a factor of ^riao ■ 

All groups with a ^ will remain 'mav — 1 Laughlin 
state with wave function (p'^ [z] . But the group with a — 
tTo will form the quasihole state 



N/s 



(21) 



The resulting state \s & vp — n fermionic liquid in 
which one fermion with spin cto has been removed at po- 
sition 77. As each boson is made out of n-|- 1 fermions, the 
corresponding bosonic state will he & u = n / (l-\-n) liquid 
in which a fraction 1/(1 -f ti) of an atom is missing in the 
component with spin do- As it follows from the theory 
of fractional quantum Hall effect, these atomic fractional 
quasiholes are anyons with 1/(1 + n) fractional statistics. 

In the analysis above we have considered a situation 
in which the laser only couples to atoms in a particular 
internal state cto. But we can also imagine situations in 
which the laser coupling affects different internal states. 
By selecting the number of internal states coupled to the 
laser we can create anyons with any fractional statistics 
of the form p/(l -l- n), with 1 < p < n. 
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The fractional statistics of the anyons we have de- 
scribed may be directly proved in an experiment similar 
to the one we proposed for detecting ^-anyons Sup- 
pose that we want to detect i-anyons. Then we will first 
prepare a system of bosons with 2 internal states, up and 
down, in a = 2/3 state. We focus then a laser, that 
only couples to atoms with spin up (or down), at position 
r/i , and increase its intensity until a single i-quasihole is 
created. Keeping constant the intensity of this laser we 
then adiabatically insert another laser at position 772, far 
enough from rji , affecting also only atoms with spin up (or 
down). The system then evolves from the one-quasihole 
state, "^ni, to a two-quasihole state, ^r]i,r]2- The crucial 
point is that at a certain point of the evolution the sys- 
tem reaches a superposition state ^' ~ 4"^^ -I- '^■rji,r]2 ■ This 
is the superposition we need to test the statistical angle. 
If we adiabatically move the laser at position 771 along a 
path enclosing position 772, the evolved state at the end 
of the process will be 

+e'2./3^^^_^^^ (22) 

where the relative phase, e*^'^/'^, reflects the statistical 
phase. This phase may be detected in a Ramsey-type 
interferometer in which the evolution of the initial super- 
position state, 5", is split in three different ways, in which 
the laser at rji performs respectively, one, two, and three 
loops enclosing 772 • If after performing the loops we con- 
tinue the evolution of the system, one can check that the 
system will only evolve to the two-quasihole state, ^r;i,772 j 
in the case in which 3 loops were performed. This 3 re- 
flects the i -statistics of the quasiparticles. 

IV. THE FERMIONIC BOSONS 

In this section we give a special attention to the case in 
which the degeneracy of spin equals the number of atoms 
in the system. Such a situation is a realistic possibility 
with a system of small number of atoms. In this case the 
effective fermions will occupy the single particle states 
with m — and spin a — 1, . . . , N . The correspond- 
ing Slater determinant is a pure spinor function, with no 
spatial dependence 

m=x[u]= en....,^«<---<, (23) 

SO that the bosonic ground state 

m] ^ x[u]Yl{z, ~ z,) (24) 

i<j 

is the product of a completely antisymmetric spin state 
and a v ^ 1 Laughlin state. We note that in the above 
wave function the spin and orbital degrees of freedom 
are completely decoupled . The orbital wave function 
is antisymmetric, so that as long as the spin degrees of 
freedom are not excited the behaviour of the system will 



be fermionic. Any property concerning orbital degrees 
of freedom that were experimentally measured would be- 
have as if we would have fermions in the system, even 
though our system is composed of bosons. 

Suppose, for example, that wc insert a laser localized 
at a certain position 77, and that the laser couples to all 
atoms, no matter the internal state. As the intensity of 
the laser is increased the system will evolve from state 
( p^ ) to the quasihole state: 

*r,[e]=xMn(^,-77)n(^,-z,)- (25) 

i i<j 

This quasihole state is a = 1 state in which one particle 
has been removed at position 77, therefore, the relation 
between filling factors and fractional statistics states that 
we have created an anyon with statistics 1; that is, a 
fermion. 



V. NUMERICAL SIMULATIONS 

In this section we illustrate the results presented in the 
previous sections with some numerical exact diagonaliza- 
tions for = 4, 6, and 8 bosons. In our calculations we 
restrict ourselves to single-particle states in the lowest 
Landau level. This approximation is valid in the limit in 
which the energy separation between Landau levels, 2a;, 
is much larger than all the energies available. 

We note that Hamiltonian (H) is invariant under rota- 
tions around the z axis, and also under global spin rota- 
tions. These symmetries allow us to diagonalize it within 
subspaces of fixed z component of the total angular mo- 
mentum, L, fixed total spin, 5, and fixed z component 
of the spin, Sz- 

Figure |] shows the eigenvalues of Hamiltonian (||) for 
a system of A^ = 6 bosons with two internal levels, in 
the rapid rotation limit. We see that there is a branch 
of states well separated by a gap from the rest of the 
spectrum. As long as the temperature is smaller than 
this gap, the system will remain within this subspace. All 
the states in this subspace lie within the lowest Landau 
level and have zero interaction energy. They form what 
we have called subspace JF. 

In order to check the fermionization scheme we have 
diagonalized the effective Hamiltonian (|^) for a system 
of A^ = 6 fermions with spin s = 1/2. Figure ^ shows 
that the bosonic spectrum projected to the subspace T is 
identical to that of the free fermions, except for a shift in 
the angular momentum of the eigenstates, AL = N{N — 
l)/2. This is precisely the angular momentum of the v = 
1 Laughlin state, that connects fermionic and bosonic 
states. For a given angular momentum, there are many 
degenerate states within the subspace JF. This density 
of states is exactly reproduced by the system of non- 
interacting fermions (Figure ||) , exept again, for the shift 
in angular momentum. 

We have studied the ground state properties of a sys- 
tem of A'^ = 4, 6, 8 bosons, and n internal states, for the 
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FIG. 1: Eigenvalues of Hamiltonian(^) for A'^ = 6 bosons and 
n = 2 internal levels, as a function of the total angular mo- 
mentum L. For illustration purposes we chose (1 — fl/uj)/g = 
0.01. The energy is measured in units of gfkJ, and all energies 
are shifted by AE — —2.1ghuj. For total angular momenta 
larger than L = 23 the dimension of the subspace is too large 
and we can not diagonalize the Hamiltonian exactly. We have 
used a variational approach to obtain the ground state and 
the first excited states. 
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FIG. 3: Density of states for A'^ = 6 bosons with n — 2 internal 
levels in the subspace J^, as a function of the total angular 
momentum L. We have also plotted the density of states for 
a system of A*' = 6 fermions and n — 2 internal levels with 
Hamiltonian (H). 



Laughlin state: 



Lo = N{N/n - 1) + N{N - l)/2. 



(26) 



We have checked that the effective fermions form in- 
teger quantum Hall states at filling factor v — n. The 
density profile of these states is nearly flat in the bulk, 
with n fermions per unit of area (Figure |). The cor- 
responding bosonic states are fractional quantum Hall 
liquids at filling factor v = n/{n + 1). The density pro- 
file of these states is also nearly fiat in the bulk, with 
n/ (1 -I- n) bosons per unit area (Figure ^). 

We have also checked that, starting from these frac- 
tional liquids, anyons with p/{l + n) statistics may be 
created. We consider first the following Hamiltonian, 



(27) 



FIG. 2: Eigenvalues of Hamiltonian (g) for TV = 6 bosons pro- 
jected to the subspace T (right), and eigenvalues of Hamilto- 
nian ^ for A'' = 6 fermions (left), as a function of the total 
angular momentum L. Both spectra are identical except for 
a shift in the angular momentum AL = A'^(A'^ — l)/2. The 
data shown correspond to (1 — Q,/ijj)/g = 0.01. The energy 
is measured in units of ghto, and all energies are shifted by 
AE = -2.1ghuj. 



cases in which N/n is an integer number. We find (Fig- 
ure ^) that the total angular momentum Lq of the ground 
state is always equal to the one of the Fermi sea formed 
by the fermions, plus the angular momentum of the v = 1 



which includes the presence of a laser localized at po- 
sition 771, affecting p different spin states. We have di- 
agonalized Hamiltonian (|2^) within the subspace for 
iV = 6 bosons, n = 1, 2, 3, 6, and 1 < p < n. When the 
laser power is sufficiently large we find that the ground 
state is a one-quasihole state, '1',,^ , in which a certain 
amount of atoms is missing in the components with spin 
states (Ti, . . . , (Tp. We have also diagonalized a Hamilto- 
nian describing the system in the presence of two lasers 
localized at positions 771 and r/2, coupled to the same p 
internal states. The resulting ground state, 4'r)i,?;2: ^^OH" 
sists of two identical quasiholes located at positions iji 
and r]2. Having the states '^rji and '^rii,ri2 '^s can check 
very easily what is the statistics of the quasiholes we have 
created. Suppose that we adiabatically drive the laser at 
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FIG. 4: Total angular momentum Lo of the ground state for 
a system with A'^ = 4, 6, 8 bosons, as a function of the number 
of internal levels n. We have plotted the numerical values 
(stars) together with the curve given by equation (Eq). 
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FIG. 5: Density profile of the ground state of a system of 
= 6 fermions with Hamiltonian (^. The different curves 
correspond to different values of the number of internal states 
n. The unit of density is 1/(2'k£^), and the unit of length is 
i. 



position 771 along a path enclosing position 772- At the 
end of the process the states vE'^j and '^rii,ri2 pick up 
Berry phases given, respectively, by. 



71 = 27r / dxdyl^^J^ 

J A 

72 = 27r / dxdy\-^,^^ 

J A 



FIG. 6: Density profile of the ground state of a system of 
= 6 bosons with Hamiltonian (P), in the limit of rapid 
rotation. The different curves correspond to different values 
of the number of internal states n. The unit of density is 
1/{2tv£'^), and the unit of length is £. 



where the integrals are performed over the area A en- 
closed by the path described by the laser. The equations 
( ^ ) for the Berry phases can be easily derived, taking 
into account that the states and are made out 

of fractional Laughlin quasihole states iQ. The differ- 
ence between the Berry phases 72 — 71 will reflect the 
extra phase that the quasihole at 771 picks up because of 
the presence of the other quasihole at 772 . Since the closed 
loop we have performed is equivalent to two consecutive 
interchanges of the quasiholes we have that the statistics 
of the quasiparticles is given by the angle 



72 - 71 



(29) 



m I 1 



(28) 



We have calculated expression (|29| ) for quasihole states 
corresponding to iV = 6 bosons and different numbers 
of internal states n = 1,2,3,6. The statistical phases 
we obtain (Figure |^) are in excellent agreement with the 
fractions p/(l + n), even though the finite size of the 
system. 



VI. EXPERIMENTAL CONDITIONS 

In this section we discuss the set of conditions that a 
system of bosonic atoms must fulfill in order to observe 
the fractional quantum Hall liquids and the anyons we 
have described in the preceeding sections. 

First of all we need a system with degenerate internal 
levels, such that the interaction between the atoms does 
not depend on the internal state. We have also made a 
two-dimensional approximation, so that the atoms must 
be confined in a trap with £ <^ £z- For the ground state 
to be a fractional quantum Hall liquid the system must 
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Finally , the most restrictive condition to observe the 
fractional quantum Hall states and their anyons is the 
temperature. In order to freeze out the excitations we 
need kT/fiio <C (1 — H/lo). For a system of = 6 par- 
ticles the observation of the states v — 1/2,3/4,3/4 re- 
quires temperatures kT/fiio <C 0.20,0.38,0.62. 



VII. THE PAIRED SUPERCONDUCTING 
STATES 
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FIG. 7: Statistical angle of the quasiholes created in a system 
of A'^ = 6 bosons, as a function of the number of internal 
states. We have plotted the numerical results (stars) together 
with the ones predicted for an infinite system (dots). The 
different points for each number of internal states correspond 
to different values of the number of internal states the laser 
is coupled to. 



be in the rapid rotation regime. From the numerical di- 
agonalizations we can calculate the critical frequency of 
rotation at which the system will be driven to a cor- 
related quantum Hall liquid. For a system of iV = 6 
particles we find that the ground state of the system be- 
comes a V = 1/2,2/3,3/4 liquid for a frequency of ro- 
tation ^l/uj > 0.80,0.62,0.38, respectively. These data 
show how the critical frequency of rotation decreases as 
the number of internal states is increased, so that the ob- 
servation of the fractional states v = n/{l + n) requires 
less faster rotating traps, as we approximate i> — 1. For 
an arbitrary number of particles N we can estimate what 
is the condition that the frequency of rotation must ful- 
fill in order for the a state v — n/{l + n) to be the 
ground state. The state — n/{l + n) has a total 
angular momentum L ^ i/~^N{N — 1) (decreasing as 
— > 1), and thus an angular momentum energy per par- 
ticle bl S — ^/uj){N — 1). In order to be in the 
rapid rotation regime we need to be much smaller than 
both the trap energy and the typical interaction energy, 
so that the conditions {l — rt/uj){N—l) ^ z/, vg/Air must 
be fulfiU. 

For creating the anyons the conditions required are 
much more demanding. We have first to focus the lasers 
within a distance ~ t. For a localization length t 1/im 
this implies an upper limit for the trap frequency of 
~ IkHz. As well, the frequency of rotation needed to 
create the anyonic excitations is larger that the critical 
frequency to observe the fractional quantum Hall liquids, 
since the total angular momentum of the quasihole states 
is 2N units of angular momentum larger. For a system 
with N — Q particles anyons with 1/2, 1/3, 1/4-statistics 
may be created for Jl/w > 0.91,0.81,0.68, respectively. 



In this section we study a situation in which the ground 
state of the atomic system is a paired state of effective 
fermions. We consider a system with and even number 
of atoms and n = 2 internal levels, up {]), and down 
(I), and we assume that the interaction between atoms 
in the same spin state is much larger than the interaction 
between atoms with different spin states. In this case the 
interaction term in Hamiltonian (^ can be approximated 
by 



N 

i<j 



pTT 



ij 



(30) 



where the operator P^j"^ projects the spin state of the 
pair of particles ij on the state aa. 

We will show that, in the limit of rapid rotation, the 
ground state of the system is the product of a, v — 1 
Laughlin state and a completely antisymmetric wave 
function. This antisymmetric wave function has the 
structure of a BCS state of fermions with p-wave pairing. 

From the results of section II it follows that in the rapid 
rotation limit the ground state, of the system must be 
both an analytic function of the z^'s and an eigenstate 
of the interaction (30) with eigenvalue 0. Let's choose 
any pair of particles i and j. For ^' to be annihilated 
by the interaction (30) we have two possibilities. Either 



the particles i and j are in a different spin state, or the 
relative angular momentum of the pair is larger than 0. 
It follows that either the spin triplet (tiij + iilj) or the 
factor {zi — Zj) is a factor of ^. How many factors of 



the form (zi 



and how many spin triplets we will 



have in the ground state will be the choice that costs the 
minimal amount of angular momentum energy, thus, we 
will have as many triplet states as possible. Note that the 
maximum number of pairs that we can arrange in triplets 
is precisely N/2, corresponding to a certain pairing of the 
particles. It follows that the ground state has the form 



^ = Pf 



i<3 



where the Pfaffian Pf is defined by 



34 ■ 



N~1,N), 



(31) 



(32) 



where Mij are the elements of an antisymmetric matrix 
and A denotes the operation of antisymmetrization, nor- 
malized such that each distinct term appears one with 
coefficient 1. 
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The state ( ^l| ) is the product of a v = 1 Laughhn 
state, in which the factor (zj — zj) appears for aU pairs 
of particles ij, and a Pfaffian factor. The Pfaffian factor 
acts in the foUowing way. It chooses a pairing of the 
particles, and for all pairs ij within the pairing it removes 
the factor {zi — zj) and substitutes it by a triplet state 
(Tiij + iitj)- III this way the system minimizes its total 
angular momentum, remaining in the subspace of zero 
interaction eigenstates. 

Since the Pfaffian factor is completely antisymmetric 
we can interpret it as the wave function of a system of N 
effective fermions, in the same way that we did in section 
II. In this case, however, the fermionic wave function is 
not an analytic function of the z^'s, and, therefore it does 
not lie within the lowest Landau level. Such a Pfaffian 
state has been shown to be the ground state solution of 
a BCS theory for fermions with p-wave triplet pairing 
p3[ . With this result we can say that the state (|l]) is 
a superconducting state of effective fermions. It would 
very interesting to derive an effective Hamiltonian for the 
fermionic wave function, in a similar way as we did in 
section II. This effective Hamiltonian should describe in 
this case fermions with attractive interactions between 
particles with different spin states. Derivation of this 
Hamiltonian will be discussed elsewhere. 



VIII. CONCLUSIONS 

In conclusion, we have shown how by rapidly rotating 
a gas of ultra cold bosonic atoms, we can drive the system 
into a strong interacting regime in which novel correlated 
ground states and excitations appear. 

We have developed a theory that exactly maps the 
problem of interacting bosons in the rapid rotation 
regime to a problem of non-interacting fermions. This 
fermionization scheme allows us to find the exact ground 



state and low-lying excitations of the system. For the 
case of a system with no internal levels we recover the 
results of our previous work, with a ^-Laughhn liquid 
as the ground state, and 1/2 anyonic excitations. In the 
presence of n internal levels and when the ratio N/ n is an 
integer number, the ground state of the effective fermions 
is a h' — n integer quantum Hall state. The correspond- 
ing bosonic ground state is a fractional multicomponent 
quantum Hall liquid with filling factor j/ = n/(l + n). 
This fractional liquid is a highly correlated state made 
out of n copies of a 1/(1 + ri)-Laughlin liquid, one for 
each internal state. Starting from these liquids we have 
shown how to create p/{l + n)-anyons, by focusing lasers 
at the desired positions, coupled top different spin states. 
The fractional statistics of these anyons may be directly 
tested in a Ramsey-type interferometer similar to the one 
we proposed to detect 1/2 anyons. 

The experimental observation of the fractional quan- 
tum Hall liquids and the anyons described in this work 
requires small temperatures, small number of atoms, and 
very fast rotating traps. The conditions required are re- 
laxed as the number of internal levels increases. A possi- 
ble scenario for the observation of these entangled states 
is a two-dimensional rotating optical lattice. As already 
shown in the recent experiment with a one-dimensional 
optical lattice it is possible to reach a regime in which a 
small fixed number of atoms is confined in each well of 
the lattice. Furthermore, since one has many identical 
copies of the system, the experimental signals are highly 
magnified. 
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